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ABSTRACT 

The large number of close-in Jupiter-size exoplanets prompts the question whether star-planet in¬ 
teraction (SPI) effects can be detected. We focused our attention on the system HD 17156, having 
a Jupiter-mass planet in a very eccentric orbit. Here we present results of the XMM-Newton obser¬ 
vations and of a five months coordinated optical campaign with the HARPS-N spectrograplR We 
observed HD 17156 with XMM-Newton when the planet was approaching the apoastron and tnen at 
the following periastron passage, quasi simultaneously with HARPS-N. We obtained a clear (« 5.5cr) 
X-ray detection only at the periastron visit, accompanied by a significant increase of the Rhk chromo¬ 
spheric index. We discuss two possible scenarios for the activity enhancement: magnetic reconnection 
and flaring or accretion onto the star of material tidally stripped from the planet. In any case, this is 
possibly the first evidence of a magnetic SPI effect caught in action. 

Subject headings: stars: individual (HD 17156) — stars: activity — stars: coronae — stars: late-type 
— X-rays: stars 


1. INTRODUCTION 

The discovery of new extrasolar planetary systems is 
today routinely carried out with a wealth of space-based 
and ground-based observational facilities, and our knowl¬ 
edge frontier has already moved on from their mere 
search to the determination of stellar and planetary prop¬ 
erties which are relevant to understand the evolution of 
these bodies during stellar lifetime. In particular, one of 
the open issues is the extent of magnetic star-planet in¬ 
teraction (SPI) and its detectable effects in systems with 
close-in massive planets (so called “hot Jupiters”). This 
phenomenon is astrophysically important for the charac¬ 
terization of planetary magnetospheres, and as a mecha¬ 
nism of energy and angular mo mentum transfer b etween 
the host star and its planet (s) ( |Cuntz et al.||2000| ). 

Observations of variable chroinospheric ernission signa¬ 
tures, such as the Ca H H&K lines cores, were reported 
for a few stars, with a mo dulation period close to the 
orbital period o f the planet (Shkolnik et al.||2005 Shkol- 
nik et al.||2008|), and interpreted as eviderice ot SPI ot 
magnetic origin. In fact, a tidally induced effect can be 
excluded because - for low-eccentricity orbits - the pe¬ 
riodic modulation is expected twice per orbital period, 
due to the two tidal bulges on the opposite sides of the 
stellar surface. Conversely, a purely stellar activity effect 
is unlikely because the relevant time scale should be the 
stellar rotation period rather than the planetary orbital 
period. 

However, significant chromospheric variability is visi¬ 
ble only at some epochs, and the existence of detectable 
SPI s ignatures at X-ray wa velengths is a matter of de¬ 
bate. Kashyap et al. (2008) showed that stars with hot 
Jupiters are statistically brighter in X- rays than stars 
without hot Jupiters, and Scharf (2010) claimed a cor¬ 
relation between the X-ray luminosity of the host star 
and the planetary mass which he ascribed to the same 
phenomenon . De taile d statistical studies by Poppen- 


( I2011P , 
ilts can 


haeger et al. (2010) and Poppenhaeger & Schmitt 
reached instead the conclusion that the above resii' 
be affected by observational biases, and that SPI induces 
quite small effects in most cases, which become impor¬ 
tant and measurable only in peculiar systems. 


^ Based on observations collected at the Italian Telescopic 
Nazionale Galileo (TNG), operated on the island of La Palma by 
the Fundancion Galileo Galilei of the INAF (Istituto Nazionale 
di Astrofisica), in the frame of the programme Global Architec¬ 
ture of Planetary Systems (GAPS). 


Table 1 
Observation Log 


Date 

(2014) 

MJD“ 

Orbital 

Phase 

^exp 

[s] 


Sep 5 

Sep 20** 

XMM-NewtonOhservations 
56905.02595 0.267-0.285 

56920.96779 0.018-0.037 

32,100 

32,900 

logLx/Tboi 
< -7.3 
-7.0 


HARPS-N Observations 



Aug 5 

56874.23161 

0.815 

900 

-5.079 ± 0.006 

Aug 7 

56876.23117 

0.909 

900 

-5.083 ±0.007 

Aug 8 

56877.21796 

0.956 

900 

-5.097 ± 0.005 

Aug 9 

56878.23348 

0.004 

900 

-5.096 ±0.003 

Aug 21^= 

56890.24703 

0.570 

900 

-5.080 ± 0.004 

Sep 9 

56909.07199 

0.457 

10,800 

-5.094 ± 0.001 

Sep 21*' 

56921.20339 

0.029 

7,200 

-5.077 ± 0.002 

Nov 10 

56971.08149 

0.380 

7,200 

-5.096 ±0.002 

Nov 25 

56986.07491 

0.087 

900 

-5.098 ±0.004 

Dec 7 

56998.03466 

0.650 

900 

-5.095 ±0.008 

Dec 14 

57005.08871 

0.983 

7,200 

-5.100 ±0.004 


Note. — a\ = ffJD — 2400000.5. b: simultaneous observations. 
High umidity night ('^ 80%). 


Rec e ntly, 


Pillitteri et al.| (2010) and Poppenhaeger &j 
) pinpointed the case of two wide binaries 


Wolk (2014) pinpointed 


where the planet-hosting primary component is much 
brighter in X-rays than expected by considering the ac¬ 
tivity level of the coeval companion. This finding sug¬ 
gests that SPI, by acting through transfer of angular mo¬ 
mentum from the planet to the star and increasing thus 
the stellar rotation and activity, may lead to enhanced 
X-ray emission and thus “rejuvenating” the star. 

Detailed studies of individual systems were also per¬ 
formed, leading to different results: the occurrence of 
repeated flaring events in HD 189733 A just after egress 
from planet occultation, observed with XMM-Newton 
(Pillitteri et al. 2011 2014 2015), suggests a system- 
atic SPI that leads to X-ray and 1'JV variability phased 
with the planetary motion; instead, a long XMM-Newton 


moni toring program of HD 179949 (Scandariato et al. 


2013), showed evidence of chromospheric and coronal 


variability most likely due to stellar rotation and to in¬ 
trinsic short-term activity evolution, but no clear signa¬ 
ture related to the orbital motion of the planet. 

From a theoretical point of view, analytical models 
were developed to explain the periodic chromospheric 
enhancemen ts as due to flare-like magnetic reconnection 
events (e.g., Lanza 2009), with an expected dependence 
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on the relative velocity between the coronal and plan¬ 
etary magnetic fields. Another key parameter is given 
by the relative strength of planetary and stellar mag¬ 
netic fields, which determines the level of intersectio n of 
the respective A lfven ic surfaces, a s investigated by |Co-| 
hen et al. (2009 1 and Lanza Sophisticated MHD 


simulations by Cohen et^ (2009) confirm that the in¬ 
teraction of stellar and planetary magnetospheres lead 
to magnetic stresses that may result in reconnection and 
reconfiguration events where energy is eventually dissi¬ 
pated by high-energy radiation. 

As a possible test of this scenario we focused our atten¬ 
tion on HD 17156, a system with a Jupiter-mass planet 
in an eccentric orbit. In fact, SPI is expected to have a 
strong dep endence on the star-planet separation (Cuntz 
.|2000| , and in systems with elongated orbits it must 


et al. 


occur preferentially or exclusively in proximity of perias- 
tron, when the stellar and planetary magnetospheres are 
close enough to interact. 

2. TARGET AND OBSERVATIONS 

The system is composed by a GO primary [B — V = 
0.64) at 75 pc from the Sun and a transiting planet 
(Mp ^ 3.2Mj) on a 21.2 days period orbit with ec¬ 
centricity e = 0.68, and periastron and apoastron dis¬ 
tances of 7.4 R* and 39.1 R^,, respec tively (semimajor 
axis a = 0.16AU, Barbieri et al.|[2009 ). The host star is 
characterized by a chromospheri c activity in dex log i?HK 
ranging from —5.06 to —5.01 (Pace 20131, typical of 
middle-aged solar-type stars. 

3. COORDINATED X-RAY AND OPTICAL CAMPAIGN 

We observed HD 17156 with the EPIC camera on board 
XMM-Newton in two visits of ^ 30 ks each, on Septem¬ 
ber 5 and September 20 2014. The aim was to observe 
the star with the planet at the periastron passage and 
far from periastron in order to seek phase-related X-ray 
variability which could be ascribed to SPI. 

The Observation Data Files (ODF) were reduced with 
SAS ver. 13.5 and the latest set of calibration hies 
(CCF). For each exposure, we obtained tables of the 
events calibrated in detector positions, arrival time, en¬ 
ergy, event PATTERN and quality flag. We Altered 
the events for the recommended flags and patterns 
(FLAG=0, PATTERN < 12), with energies in the 
broad band 0.3-8.0keV, and in a soft band 0.3-1.5keV. 
Since a soft spectrum is expected for low-activity stars 
like HD 17156, we preferred the soft band for source de¬ 
tection in the EPIC images, as already experimented 
with the observations of HD 1 89733, having the bu lk of 
X-ray emission below 1.5 keV (Pillitteri et al. 2014). 

For the detection of the star and for estimating signifi- 
cance and rate, or upper limits, we used a code based on a 
multi-scale wavelet convolution (|Damiani et al.|1997a|b|). 

HD 171 56 was also observed w ith the HARPS-IN spec- 
trograph (Cosentino et al.|2012 ) at the Italian 3.6m Tele¬ 
scopic Nazionale Galileo. In this paper we present data 
acquired between 2014 August and December, in prepa¬ 
ration and as a follow-up of the XMM-Newton program. 
The full observation log is shown in Tab.ffl Note in 
particular two HARPS-N observations coordinated with 
XMM-Newton, on 2014 September 9 and 21, the lat¬ 
ter effectively simultaneous with XMM-Newton while the 
planet was at periastron. 


Optical spectra were acquired with different exposure 
times, for different aims: short exposures (15 minutes 
each) belong to a long-term series aimed to precise radial 
velocity measurements, while long exposures (120 or 180 
minutes) were optimized for high-resolution spectroscopy 
of the Ca H H&K 3968.5, 3933.7 A chromospheric emis¬ 
sion lines. In order to avoid saturation at longer wave¬ 
lengths, the long exposures were always split in a number 
of 15 minutes sub-exposures, and the raw images were 
combined a posteriori to compute an average spectrum 
using the standard HARPS-N pipeline. In this respect, 
we have followed the same observa tion strategy and d ata 
reduction and analysis detailed in Borsa et al. (2015). 

These data were employed to update the orbit m pa¬ 
rameters. Adopti ng the transit timing reported in |Nutz-| 
man et al.| (|2011 ), we obtained the following ephemeris 


(RJDtdb) for the passage at periastron: 


To = (2454884.3106±0.0050)-f(21.216398±1.6xl0-®)T. 

( 1 ) 

The orbital phases at the start of each observation are 
reported in Tabj^ 


4. RESULTS 

4.1. X-Ray Emission Switching On 

In Fig. nwe show the field of view around HD 17156 in 
the two XMM-Newton exposures. The brightest source 
visible in both images is an unidentified background 
source, at ~ 10.5” from the optical position of our tar¬ 
get. HD 17156 was not detected while the planet was 
approaching the apoastron (September 5), at a level of 
3cr of the local background, while it was detected in 
the visibility window that started 9 hr after periastron 
(September 20) and lasted for about 10 hr. Employing 
the PN and MOS 2 datcQ we obtained a 5.5cr significance 
level in the soft band, with a PN-equivalent count rate of 
(7.1 ± 2.5) X 10“"^ cts s“l using a maximum wavelet spa¬ 
tial scale of 5.6”. About 26 ± 9 photons were collected in 
total in the two EPIG detectors, and their time distribu¬ 
tion is compatible with a constant source (Kolmogorov- 
Smirnov test probability > 50%). 

The evaluation of the target X-ray flux is hampered by 
the presence of the contaminating background source, de¬ 
tected in the first observation with a count rate ~ 50% 
higher than in the second observation. However, this 
source is quite hard and highly absorbed: in fact, its 
spectrum can be described as a thermal source with 
kT ^ 1.6 keV attenuated by a hydrogen column density 
Nh = 1.3 X 10^^ cm“^ . Hence, we employed only the 
events collected with the PN detector in the soft band 
(0.3-1.5 keV) already used for the source detection pro¬ 
cess. Taking into account the Encircled Energy Frac¬ 
tion for a point source within a 5.6” radius, we obtain a 
count rate of (4.8 ± 2.5) x 10“"^ cts/s for HD 17156. We 
estimated the source X-ray flux assuming optically thin 
emission from a thermal plasma with solar composition 
and T = 2 MK, as appropriate for a solar-type corona, 
thus obtaining /x = 5.7 x 10“^® erg cm“^ s“^ . A similar 
procedure yielded an upper limit /x < 3.3 x 10“^® erg 
cm“^ s“^ for the source emission far from the periastron. 

^ MOS 1 data were discarded because of the presence of a CCD 
gap near the target location. 
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Figure 1. X-ray images of HD 17156 taken far from the planetary periastron (left side panels) and near the periastron (right side panels). 
Top row panels are intensity images, bottom row panels are RGB images with photon energy colors R = 0.3 — 1.0 keV, G = 1.0 — 2.5 keV, 
B = 2.5 — 5.0 keV. Smoothing is applied to the images, with a Gaussian of tr = 2.4" for the intensity images and a = 4" for the RGB 
images. Positions of the only two objects in the SIMBAD catalog are shown with small squares. Gircle sizes indicate the wavelet detection 
scales of HD 17156 and of an unrelated background object with an harder spectrum. 


-5.06 



-5.1 1 


-5.17 1 I . 1, , I. . ,1 

0 50 100 0.0 0.2 0.4 0.6 0.8 1.0 

Day [MJD - 56874] Phase 

Figure 2. Ca H H&K chromospheric emission index with 3cr 
error bars vs. time (left) and vs. phase (right). Symbol sizes indi¬ 
cate different exposure lengths. In red the HARPS-N observations 
coordinated with XMM-Newton, and the baseline obtained with 
an error-weighted regression, excluding only the data point which 
corresponds to the X-ray detection of the target. 

In order to evaluate the intrinsic source X-ray lumi¬ 
nosity, we estimated the interstellar absorption starting 



we obtain an unabsorbed flux of /x = 1.5 x 10“^^ 
erg cm“^ s“^, and a corresponding X-ray luminosity of 
Lx ^ 1.0 X 10^^ erg s“^ close to periastron. Similarly, 
far from periastron the upper limit is < 6.0 x 10 ^® erg 
s“ ^ . We note however that these estimates are very sen¬ 
sitive to the assumed plasma temperature. For example, 
a 1 MK plasma would imply X-ray luminosities 6 times 
larger than reported above. 


4.2. Chromospheric emission 

The Mount Wilson i?HK index was comp uted by the 
HARP S-N pipeline following the recipe of Lovis et al. 
(2011|). We recall that this is a measure of the chro¬ 
mospheric emission in the Ca II H&K lines, corrected 
for a basal level depending on the stellar B — V color, 
and normalized to the stellar bolometric luminosity. In 
Fig. I^we show the variation of logi?^^: time, from 

2014Wugust to December, and vs. orbital phase. In the 
5 months considered, the data suggest a chromospheric 
activity level essentially flat within the uncertainties of 
the measures. Most remarkable is the significantly higher 
(> 8 cr) value of Rhk nieasured simultaneously with the 
X-ray detection of our target. In fact, we observed a 
4% increase of the chromospheric emission at phase 0.03 
(just after periastron) with respect to the previous mea¬ 
surement taken at phase 0.46 (approaching apoastron). 
High values apparently occurred also on other dates 
in 2014 August, but these measurements are affected by 
relatively large uncertainties because they were derived 
from observations with short exposures. Some low-level 
variability could also be due to resi dual systematics i n 
the reduction of the HARPS-N data ( Lovis et al.||2011 ). 


5. DISGUSSION AND GONCLUSIONS 

HD 17156 showed enhanced chromospheric and coro¬ 
nal emission a few hours after the passage of the planet 
at the periastron. Both the X-ray luminosity and the 
i?(jK index indicate a star with a magnetic activity 
level lower than expected. In fact, modeling of the 
Rossiter-McLaughlin effect yields a rot ational velocity 
of 4.2 ± 0.3km s“l (Narita et al. 2009), and the com¬ 
puted rotation period Prot ~ IS days leads to a presumed 
X-ray luminosity = 1-2 x 10^^ erg s~^, adoptin g 


the activity-rotation relations by Pizzolato et al. (2003). 
Starting from a stellar age of 3.37j'jn 47 Gyr, based on 
transit and asteroseismic observations by |Nutzman et ^ 
(2011), we also derive ~ —4.8 and Lx/Lhoi i n the 
range [—5.8, —6.0] from the activity-age relations by Ma- 
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majek & Hillenbrand (2008). Our observed chromo- 
spheric and coronal emission levels are, respectively, a 
factor ~ 2 and ~ 10 lower than the values above. On 
the other hand, the measured chromospheric index in¬ 
dicates an X-ray luminosity 1-2x10^^ erg s“^, in agree¬ 
ment with the XMM-Newtonohseiva.tio'n.s. We conclude 
that HD 17156 was in a low-activity state, with a “qui¬ 
escent” X-ray luminosity just below the detection sen¬ 
sitivity of the Sept 5 XMM-NewtonohseTvaXion. Con¬ 
sequently, in the following we assume an average large- 


scale magnetic field strength of just IG (Vidotto et al. 
2014), and an energy release AL^ ~ 5 x 10^*^ erg s“^ due 


to an SPI event at the periastron. We deem unlikely 
the occurrence of a spontaneous flaring event of such 
power (ALx ~ quiescent L^) in the corona of a quiet 
star, such as HD 17156, just after reaching the minimum 
star-planet separation. 

One hypothesis is that the SPI phenomenon has a mag¬ 
netic origin. We evaluate d the energy budget with the 
magnetic SPI models by |Lanza| (|2013D . We assume a 
loop interconnecting the stellar magnetic field with a 
planetary field of ^ 10 G (based on analogy with our 
Jupiter and theoretical models, e.g. Reiners fc Ghris- 


tensen 


2010). For magnetic confinement, the plasma 


temperature near the planet position should not exceed 
0.7 MK. In this case, the maximum available power due 
to magnetic stresses near periastron is Pstress ~ 2 x 10^^ 
erg s“^, for a stellar potential dipolar configuration. We 
do not expect any energy release by this mechanism at 
apoastron. A possible difficulty with this explanation 
rests on the MHD time scales for establishing an inter¬ 
connecting loop, building up the magnetic stresses and 
releasing the accumulated energy, given the fast plane¬ 
tary passage. In fact, the relative velocity of the mag¬ 
netic field co-rotating with the star with respect to the 
planet speed at periastron is about 160 km s“^, implying 
a dynamical time scale of 10-15 hr. Stability of such a 
long loop {I ~ 10^^ cm) is also an open issue. Recon¬ 
nection between separate stellar and planetary magne¬ 
tospheres is also unlikely, because the dissipated power 
would be lower tha n estimated above by about four or¬ 
ders of magnitude ( Lanza|[2012 ). 

An alternative possibility is that the power emitted in 
X-rays is released by matter evaporated from the planet 


and accreted onto the star nea r periastron (Matsakos 


et al. 2015 Pillitteri et al. 2015). This requires a mass 
accretion rate larger than ~ 3.6 x 10^^ g s“^ to make a 
power of 5 X 10^° erg s“^ available. The free-fall time 
from the inner Lagrangian point LI to the star is about 
14 hr, making this process compatible with the delay of 
the brightening we observed after periastron. On the 
other hand, detailed simulations are required to under¬ 
stand the mechanism and time scale of accumulation and 
spill over of evaporated material onto the star. 

In any case, inspection of Fig. shows that no signifi¬ 
cant chromospheric brightening was observed near other 
periastron passages, preceding or following the one high¬ 
lighted by the X-ray detection. We conclude that SPI 
events in this high-eccentricity system occur occasion¬ 
ally, or they have a short time duration. Monitoring of 
HD 17156 will be continued in the frame of the GAPS 


programme with the aim to better understand the fre¬ 
quency of the putative SPI effects and their nature. On 
the other hand, new X-ray observations are severely lim¬ 
ited by the faintness of the target and by visibility con¬ 
straints with XMM-Newton. 


AM and the GAPS project team acknowledge support 
from the “Progetti Premiali” funding scheme of the Ital¬ 
ian Ministry of Education, University, and Research. IP 
acknowledges support from the European Union Seventh 
Framework Programme FP7/2007-2013, under the grant 
agreement no 267251 “Astronomy Fellowships in Italy” 
(AstroFIt). 

Facilities: XMM-Newton{FPlC), TNG (HARPS-N) 
REFERENCES 


Barbieri, M., Alonso, R., Desidera, S., et al. 2009, A&A, 503, 601 
Borsa, F., Scandariato, G., Rainer, M., et al. 2015, A&A, 578, 
A64 

Cohen, O., Drake, J. J., Kashyap, V. L., et al. 2009, ApJ, 704, L85 
Cosentino, R., Lovis, C., Pepe, F., et al. 2012, in Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 8446, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, 1 
Cuntz, M., Saar, S. H., Sz Musielak, Z. E. 2000, ApJ, 533, L151 
Damiani, F., Maggio, A., Micela, G., Sz Sciortino, S. 1997a, ApJ, 
483, 350 

—. 1997b, ApJ, 483, 370 

Gilliland, R. L., McCullough, P. R., Nelan, E. P., et ah 2011, 

ApJ, 726, 2 

Giiver, T., Sz Ozel, F. 2009, MNRAS, 400, 2050 

Kashyap, V. L., Drake, J. J., & Saar, S. H. 2008, ApJ, 687, 1339 

Lanza, A. F. 2009, A&A, 505, 339 

—. 2012, A&A, 544, A23 

—. 2013, A,^A, 557, A31 

Lovis, C., Dumusque, X., Santos, N. C., et al. 2011, ArXiv 
e-prints, arXiv:1107.5325 

Mamajek, E. E., Sz Hillenbrand, L. A. 2008, ApJ, 687, 1264 
Matsakos, T., Uribe, A., Sc Kdnigl, A. 2015, ASzA, 578, A6 
Narita, N., Hirano, T., Sato, B., et al. 2009, PASJ, 61, 991 
Nutzman, P., Gilliland, R. L., McCullough, P. R., et al. 2011, 
ApJ, 726, 3 

Pace, G. 2013, A&A, 551, L8 

Pillitteri, L, Gunther, H. M., Wolk, S. J., Kashyap, V. L., Sc 
Cohen, O. 2011, ApJ, 741, L18 
Pillitteri, L, Maggio, A., Micela, G., et al. 2015, ApJ, 805, 52 
Pillitteri, L, Wolk, S. J., Cohen, O., et al. 2010, ApJ, 722, 1216 
Pillitteri, L, Wolk, S. J., Lopez-Santiago, J., et al. 2014, ApJ, 785, 
145 

Pizzolato, N., Maggio, A., Micela, G., Sciortino, S., Sz Ventura, P. 
2003, A&A, 397, 147 

Poppenhaeger, K., Robrade, J., Sc Schmitt, J. H. M. M. 2010, 
A&A, 515, A98 

Poppenhaeger, K., Sc Schmitt, J. H. M. M. 2011, ApJ, 735, 59 
Poppenhaeger, K., Sc Wolk, S. J. 2014, A&A, 565, LI 
Reiners, A., Sz Christensen, U. R. 2010, A&A, 522, A13 
Scandariato, G., Maggio, A., Lanza, A. F., et al. 2013, A&A, 552, 
A7 

Scharf, C. A. 2010, ApJ, 722, 1547 
Sekiguchi, M., Sc Fukugita, M. 2000, AJ, 120, 1072 
Shkolnik, E., Bohlender, D. A., Walker, G. A. H., Sc Collier 
Cameron, A. 2008, ApJ, 676, 628 
Shkolnik, E., Walker, G. A. H., Bohlender, D. A., Gu, P.-G., Sc 
Kiirster, M. 2005, ApJ, 622, 1075 
Vidotto, A. A., Gregory, S. G., Jardine, M., et al. 2014, MNRAS, 
441, 2361 























